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Classical optics can be used to efficiently implement certain quantum information processing tasks 
with a high degree of control, for example, one-dimensional quantum walks through the space of 
orbital angular momentum of light directed by its polarization. To explore the potential of quan¬ 
tum information processing with classical light, we here suggest a method to realize d-dimensional 
quantum walks with classical optics—an important step towards robust implementation of certain 
quantum algorithms. In this scheme, different degrees of freedom of light, such as frequency, orbital 
angular momentum, and time bins, represent different directions for the walker while the coin to 
decide which direction the walker takes is realized by employing the polarization combined with 
different light paths. 

PACS numbers: 42.50.Ex, 05.40.Fb, 42.50.Tx 


Introduction. A quantum algorithm is an ordered set 
of instructions given to a quantum computer to realize a 
particular processing of quantum information. However, 
due to the delicate nature of quantum systems, it is not 
easy to construct a large enough quantum computer that 
is capable of implementing sufficiently complex quantum 
algorithms. Interestingly, a number of quantum algo¬ 
rithms, such as the Deutsch algorithm [TJ, the Deutsch- 
Jozsa algorithm Q, and Grover’s search algorithm 0, 
can be simulated and implemented with certain classical 
systems Q among them with classical states of light [gi]. 


In principle, quantum walks can be used to realize 
quantum algorithms Algorithms such as quan¬ 

tum search Toj, U | and graph isomorphism 12 -14 have 
already been demonstrated using quantum walks. Re¬ 
cently, it was proved that the most general quantum 
algorithms can be implemented using only quantum 
walks EM, thus making them universal for quantum 
computing. However, there is no efficient way known to 
implement an arbitrary algorithm using a quantum walk 
with a single walker. 


Quantum walks can be realized with classical sys¬ 
tems 13 - 22 } that are simpler to prepare, control and 
measure than their quantum counterparts. Despite all 
these advantages, so far, only one-dimensional quantum 
walks have been implemented 21, 23-28}. There have 
been a few attempts to realize two-dimensional quantum 
walks experimentally j22j. However, since the complexity 
of the realization schemes for two-dimensional quantum 
walks is much higher than that of their one-dimensional 
counterparts, it is unlikely that we see three- or higher¬ 
dimensional quantum walks merely based on quantum 
systems. For example, in Ref. 22], Schreiber et al. 
used a time-multiplexing technique with single photons 


to implement two-dimensional quantum walks. In this 
method, a two-dimensional lattice is mapped onto the 
one-dimensional discrete time line. On the other, hand 
Roldan et al. E3 presented a scheme to implement two- 
dimensional quantum walks where the propagation of the 
walker in the x and y directions is the displacement in 
the frequency domain corresponding to their respective 
linear polarizations. 


In this Rapid Communication, we present a scalable 
optical implementation scheme for a particular class 
of multidimensional quantum walk. In this class a 
d-dimensional quantum walk is equivalent to d one¬ 
dimensional quantum walks alon g w ith a 2d-dimensional 
coin operator mm (see also ]3lJ). The scheme pre¬ 
sented here requires only a linear optical setup and is 
capable of performing d-dimensional quantum walks on 
a regular lattice, where d is the number of independent 
degrees of freedom (DoFs) of paraxial light. It can be 
implemented using a single photon but also with coher¬ 
ent classical light. This is due to the analogy between 
the optics of single photons and coherent states prepared 
in the same light modes. In particular, the ingredients 
of quantum walks, namely, superposition, interference, 
and, indeed, a form of entanglement, are also present 
in classical optics [32]]. This analogy is used to simu¬ 
late quantum walks for single photons with weak coher¬ 
ent states of light [U [22}. In Ref. m , we have argued 
that one-dimensional quantum walks with a single walker 
can be implemented using classical light. Here, we claim 
that classical light can also be used for multidimensional 
walks. Although possible, there are disadvantages to us¬ 
ing single photons to implement multidimensional quan¬ 
tum walks compared to classical light, as we discuss in 
the conclusions. 
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In general, a d-dimensional quantum walk results in 
multipartite entanglement between the different coordi¬ 
nates (as independent infinite-dimensional quantum sys¬ 
tems) and the coin. In the context of a realization with 
classical light this would amount to nonquantum entan¬ 
glement between the different DoFs of the light beams. 
Here, we present a scheme where the classical counterpart 
of multipartite entanglement can be observed. 

To start with, consider a one-dimensional quantum 
walk with a two-dimensional coin. The dynamics of the 
walker can be characterized by two operators: (i) the 
conditional shift operator S and (ii) the coin toss oper¬ 
ator U c . Let |t) and ||) be two orthogonal states of the 
coin. Then, the shift operator S for the walker on the 
line can be defined as 

5 = F®|t)(t|+F t ®|i)(t|. (1) 

Here, F = Y2j I j + 1) 01 and F t represent forward and 
backward propagators on the lattice, respectively. The 
coin toss operator U c shuffles the state of the coin after 
every shift, which can be defined by its action on the coin 
states of the walker as 

It) — > a It) + P ID > It) “* ~P* It) + a * It) ) (2) 

and it belongs to the group of 2 x 2 unitary matrices 
SU(2). Here, a and P are two complex numbers such 
that |a| 2 +1/3| 2 = 1. Thus, the quantum walk propagator 
reads 

Z = S( 1 ® U c ) =F ® |t) (a (tl - P* (tl) 

+ -fd <8> |t) (P (t| + a * (t|)- (3) 

The combined operation Z causes a toss of the coin state 
and a coin-state-dependent propagation of the walker on 
the lattice. 

One can extend the idea of one-dimensional quantum 
walks to multiple dimensions. The conditional shift op¬ 
erator [cp. Eq. ([l])] for quantum walks in d dimensions 
reads 

d 

s ( 2d ) = E (^nSltnXtnl+i^SltnHtnl), (4) 

n—1 

where F n and F£ correspond to the forward and back¬ 
ward propagators for the nth direction and the sets 
{|tn) ! 14-n)} constitute the basis for a 2d-dimensional 
coin space. The coin toss operator, on the other hand, 
is a unitary operator Uc 2d ^ £ U(2d) acting on the 2d- 
dimensional state space of the coin. Thus, the quantum 
walk propagator Z 2d for d-dimensional quantum walks 
reads 


Z 2d = S { - 2d) {t®U? d) ). (5) 


Note that the conditional shift operator S^ 2d ' ) defined 
in Eq. © can be rewritten as the sum of operators S n as 


S ( 2d ) = j2 s ^ 


( 6 ) 


n= 1 


where each of the operators S n represents a one¬ 
dimensional conditional shift operator identical to the 
one defined in Eq. ©. Thus d-dimensional quantum 
walks can be realized by using d independent, one¬ 
dimensional conditional shift operators assisted by a 2d- 
dimensional coin toss operator. Therefore, the main chal¬ 
lenges in realizing a d-dimensional quantum walk are (a) 
devising methods to implement d one-dimensional inde¬ 
pendent quantum walks and (b) methods to realize 2d- 
dimensional coin operations. Classical light beams pos¬ 
sess a number of independent DoFs such as orbital an¬ 
gular momentum (OAM), time spacing, frequency spac¬ 
ing, and spatial displacements. There are several schemes 
to implement one- and two-dimensional quantum walks 
in these DoFs [H-ED, m, thus establishing that clas¬ 


sical light satisfy the first requirement (a). In the fol¬ 
lowing, we demonstrate a scheme to realize an arbitrary 
2d-dimensional coin operation in classical light using only 
linear optics. 

Coin toss operator. In our scheme, the 2d-dimensional 
coin is realized by using d copropagating light beams, 
each possessing a 2-dimensional polarization space. In 
such a realization, the unitary operations on polariza¬ 
tion states in individual beams and the beam split¬ 
ters to mix different beams constitute the set of el¬ 
ementary operations. To construct an arbitrary 2d- 
dimensional coin operation Uc 2d ' 1 based on linear optics, 
we need to decompose it into these elementary opera¬ 
tions. The CS (“cosine-sine”) decomposition for unitary 
matrices fH, 3] enables us to achieve this. 

For example, using this decomposition we can write 
any unitary operator acting on a four-dimensional 
complex Hilbert space as 


with 


UJ (4) = 


( L \ 

° 

=1 

° 
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(7) 


( 8 ) 


The operators i?i, R 2 , Li, and L 2 are two-dimensional 
unitary matrices acting on the polarization in each beam. 
These operators can be realized by using a combina¬ 
tion of quarter-wave-plates and half-wave-plates 35, 
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(a) 


and I 2 represent the half-wave plates with the fast axis 
in the horizontal direction and a 2 -dimensional identity 
operator, respectively (see Fig. ©)• 

The decomposition of the four-dimensional coin op¬ 
erator Uc given in Eq. 0 can be extended to 2 d- 
dimensional coin operators [34|. In that case the ma¬ 
trices L 1 and R± are two-dimensional unitary operators 
whereas, L 2 and i ? 2 are (2d — 2 )-dimensional unitary ma¬ 
trices. The matrix can also be implemented by a 

polarization sensitive beam splitter acting on two beams 
and is of the form: 



(b) 


FIG. 1. (Color online) (a) Diagram for the optical imple¬ 
mentation of a four-dimensional unitary operation. Here, BS 
stands for the elementary beam-splitter array implementing 
the operator U*© in the CS decomposition 0 , which consists 
out of two beam splitters and two half-wave plates and shown 
in (b). 



-*» 

c 


FIG. 2. (Color online) Any unitary coin toss operator for a 
quantum walk in three dimensions can be implemented by 
a combination of three elementary beam-splitter arrays U © 
combined with 9 unitary transforms of the polarization. 


whereas the matrix UJ© is a four-dimensional unitary op¬ 
erator that is implemented by means of two beam split¬ 
ters and two half-wave plates, 

U (4) = (B(9 3 ) (g) 1 2 ) (1 2 © o z ) (B(0i) © 1 2 ) (1 2 © <Jz ), 

(9) 

where 

Bm = ( cose 0 do) 

\ sm 9 cos 0 I 

represents the beam-splitter matrix with transmittivity 
cos 2 9 , 9 3 = (9i+9 2 )/2, and 9 4 = (9\ — 0 2 )/2, whereas cr~ 




U© 

0 

0 

1 

-e 

<M 


( 11 ) 


The operators L 2 and R 2 can further be decomposed 
in 2-dimensional and (d — 4)-dimensional operators, and 
hence using CS decomposition repetitively we can create 
U™ by means of the elementary operations. In Fig. [2] 
we present a schematic representation of this decompo¬ 
sition for d = 3. Thus, by using only the elementary 
operations we can perform an arbitrary 2 d-dimensional 
unitary operation. 

Once we have a 2d-dimensional coin operator the 
rest of the scheme to implement the corresponding d- 
dimensional quantum walk is fairly straight forward. In 
Fig. [3] we present the schematics for the realization of a 
three-dimensional quantum walk. In this figure we show 
a coin toss operator Uc followed by three shift opera¬ 
tors Si, S 2 , and S 3 in three independent paths named 
a, b and c. The shift operators Si, S 2 and S 3 (in gen¬ 
eral) act on independent DoFs. For example, Si can be 
chosen as the shift operator acting on OAM states of 
the light beam • Likewise, the operators S 2 and S 3 
could be implemented as shifts of the walker in time or 
in frequency domains 19!, 21 > 22 | and in position space. 
Hence, using our scheme, it is possible to implement a 
three-dimensional quantum walk with classical light. 

Quantum walk with Grover’s coin. As an example, in 
this section we present the scheme for qua ntum walks 
in two dimensions using the Grover coin [371 ] . which, for 
two-dimensional quantum walks, is given by 


[A 4) = i 


/-1 1 1 i\ 

I - 111 

II- 11 
V 1 1 1 - 1 / 


( 12 ) 


Using the CS decomposition given in Eq. 0, this oper¬ 
ator can be decomposed as 


U, 
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(13) 
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where H and Ug read 


H = 



U G = 


(° 

0 

-1 

°\ 

0 

1 

0 

0 

1 

0 

0 

0 

v° 

0 

0 

1/ 


(14) 


The operator H can be realized by a half-wave plate with 
the fast axis at an angle n/8 with respect to the horizon¬ 
tal axis, whereas the operator IUg can be realized using 
the scheme given in Eq. © with 63 = O 4 = 7 t/4 . 

If we choose OAM and the time bins as two indepen¬ 
dent DoFs, then any state vector of the walker can be 
written as a generalized Jones vector (path-polarization 


vector |38H41i]). 


( E t, tn a t,t n E(t n ) exp(i£(j)) ^ 
Eqt„ Pe,t n E{t n ) exp(i£(j)) 
Er,t„ lt,t n E(t n ) exp(i£(j)) 

V E*,t„ h,t n E{t n ) exp (ity) ) 


(15) 


where a, j3 , 7 , <5 are complex amplitudes of 
horizontally and vertically polarized light in 
the first and the second beam such that 

E*, tn (Kd 2 + Ift.d 2 + b, t J 2 + IdJ 2 ) = 1, 

whereas E(t n ) denotes the temporal shape of the laser 
pulse centered around time t = t n . It can be chosen to 
be a Gaussian pulse with temporal width r and can be 
written as 


E(t n ) = exp - 


(t - t n f 


2 T 2 


In this notation the shift operators read 


Si = 


( e i(t> 0 0 0 \ 

0 0 0 
0 0 0 0 

\ 0 0 0 0 / 


,s 2 = 


( 0 0 0 
0 0 0 
0 0 /(At) 
\0 0 


(16) 


0 

0 


0 /(-At)/ 
(17) 


Here, At is the step length in the time domain and /(At) 
represents an operator such that 


f{At)E(t n ) = E(t n + At), 


(18) 


which can be realized by a Mach-Zehnder interferometer 
with unequal path lengths. 

If the initial value of the OAM of the walker is zero, 
his initial time is to and the light is initially travelling 
in the first beam with horizontal polarization, then the 
state of the walker reads 


( 1 \ 
0 
0 
V 0 



FIG. 3. (Color online) Schematic diagram for the setup for 
a 3-dimensional quantum walk. Here, a, b and c denote the 
three paths for the light beams, U c is the 6 dimensional coin 
flip operator, and Si, S 2 , and S 3 are the three shift operators 
acting on OAM, time, and the position space, respectively. 


After a single step of the quantum walk the state trans¬ 
forms to 


IVd) = su^ 1^0) 


1 

2 


/ - E{t 0 )e* \ 
E(to)e~^ 
E(to + At) 

V E{t 0 - At) 


m 


Thus, by repeating this operation, one can perform the 
Grover quantum walks in two dimensions. 

Multidegree nonquantum entanglement. The scheme 
presented in this Rapid Communication enables a mul¬ 
tidimensional quantum walk in different DoFs of light. 
The classical nature of the scheme provides the usual 
benefits such as immunity against disturbance, it is tech¬ 
nically simple to realize, and measurements can track the 
spread of the walk in real time. Beside these advantages, 
our scheme provides a platform to realize the classical 
counterpart of multipartite quantum entanglement, i.e., 
multidegree entanglement in classical light beams. This 
is due to the entangling property of quantum walk evolu¬ 
tion, which tends to entangle different coordinates of the 
walker. In our scheme the independent DoFs of the light 
beam serve as different coordinates of the walker. Thus, 
the quantum walk evolution results in entanglement be¬ 
tween these DoFs of the classical light beam. However, 
since the whole setup is classical in nature, the entangle¬ 
ment in question is not quantum j 20 i . l4fj | . 

To understand this in more detail, let us consider again 
the Grover quantum walk. After a single step of the 
quantum walk the state of the walker [cf. Eq. (l20l) ] can 
be rewritten as 

IV’i) = 2 ( I 1 ) l*o) l^i) + l —1 ) \to) K) 

+ |0) 11 0 + At) \h 2 ) + |0) 11 0 — At) ^ 2 )), (21) 


IV’o) = E(t 0 ) 


(19) 
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where h\,v\,h 2 ,V 2 represent the horizontal and vertical 
polarization inodes in the first and the second beam and 
|1) \to) \hi) = [ E(to)e ^ 0 0 0] T , etc. It is clear from 
Eg. (1211) that the state |i/>i) cannot be written as the 
product of OAM modes, time bins, and the polariza¬ 
tion, hence it represents an entangled state. However, 
since the whole setup is classical, the entanglement is 
nonquantum. 

In conclusion, we have presented a scheme to real¬ 
ize multidimensional quantum walks using the different 
DoFs of a light beam. The scheme uses only bright classi¬ 
cal light and is therefore very robust. Bright laser pulses 
are considerably simpler to produce, to process, and to 
detect than single photons. This makes it possible to 
perform a quantum walk on a large lattice. With the 
state-of-the-art technolo gy i t is possible to measure 100 
OAM states of light 2di 441. Similarly, one-dimensional 
quantum walks can be performed for 50 steps in time 
bins [2l|,[22| and for 10 steps in spatial modes [27j. Thus, 
using our scheme and these available methods to perform 
one-dimensional quantum walks, we can perform between 
10 and 20 steps three-dimensional quantum walks. More¬ 
over, in principle, losses can be compensated by ampli¬ 
fication unlike in the quantum regime of light, resulting 
in many steps possible in the multidimensional quantum 
walk. One of its most intriguing aspects is the occurrence 
of the classical counterpart of quantum multipartite en¬ 
tanglement, i.e., entanglement between more than two 
DoFs. Since quantum walks are not possible without the 
entanglement between different directions along with coin 
states, multi-degree entanglement is essential for classical 
implementation of the multidimensional quantum walks. 
Devising the use of this particular type of entanglement 
will motivate our future studies. 
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